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Some chemical reactions are described by electron transfer (ET) processes. The underlying mech- 
anism could be modeled as a polaron motion in the molecular crystal — the Holstein model. By 
taking spin degrees of freedom into consideration, we generalize the Holstein model (molecular crys- 
tal model) to microscopically describe an ET chemical reaction. In our model, the electron spins in 
the radical pair simultaneously interact with a magnetic field and their nuclear-spin environments. 
By virtue of the perturbation approach, we obtain the chemical reaction rates for different initial 
states. It is discovered that the chemical reaction rate of the triplet state demonstrates its depen- 
dence on the direction of the magnetic field while the counterpart of the singlet state does not. This 
difference is attributed to the explicit dependence of the triplet state on the direction when the axis 
is rotated. Our model may provide a possible candidate for the microscopic origin of avian compass. 

PACS numbers: 34.70.+e, 03.65.Xp, 82.30.-b, 82.39.Jn 



I. INTRODUCTION 

Nowadays, it has been prevailing in both experimen- 
tal and theoretical explorations that quantum coherence 
effect due to the role of phase in quantum superposition 
may exist in living processes. This essentially implies 
that there may exist quantum coherence effect in chemi- 
cal reactions in some living processes, such as charge and 
energy transfer in photosynthesis (j-jg and singlet-and- 
triplet transition in avian compass (6l4l6|. 

It has long been questioned how migratory birds can 
navigate to their destination over hundreds of miles. One 
of the possible answers is given by the radical pair mech- 
anism |9l-[lll|. Two unpaired electron spins in the radical 
pair are initially prepared in the singlet state. Due to 
their interactions with the geomagnetic field and their en- 
vironmental nuclear spins, the election spins coherently 
transit between the singlet and triplet states. Since the 
singlet and triplet states could result in different prod- 
ucts of chemical reactions, the direction and magnitude 
of the geomagnetic field determine the relative yields of 
two distinct products. By sensing the information incor- 
porated in the final products of the chemical reactions in 
their retinas, the birds can find their way to their des- 
tination. Therefore, the quantum coherence underlies in 
the avian compass since the singlet and triplet spin states 
correspond to different types of quantum entanglement. 
Ever since it was proposed a decade ago, the radical-pair- 
mechanism-based avian compass has been in favor by a 
series of biological and chemical experiments [6, 7]. 

In this hypothesis, the nuclear spins play a crucial role 
because there would be no coherent transition between 
the singlet and the triplet states if there were no nuclear 
spins [ll|. Previous studies mainly concentrated on the 
nuclear-spin environment without inter-coupling fl3l-[l5j. 
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Mostly recently, by taking into account the inter-coupling 
of the nuclear spins, we studied a special avian compass 
model with the nuclear environments modeled by an Ising 
model in a transverse field [16]. The rationality of this 
model lies in the fact that the weak inter-nuclear-spin 
coupling is comparable with the Zeeman energy splitting 
induced by the weal geomagnetic field. It was discov- 
ered that the quantum criticality in the environments en- 
hances the sensitivity of magneto-reception. On the other 
hand, although various master-equation approaches were 
proposed to deal with such spin-dependent chemical re- 
actions in the avian compass [7], the underlying physi- 
cal mechanism is still missing in studying the quantum 
coherence with microscopic models. Thus, it is urgent 
to propose appropriate microscopic models for different 
kinds of chemical reactions to make the quantum coher- 
ence effect in those processes better understood. A case 
in point is the Holstein's molecular crystal model, which 
is also regarded as a microscopic model of chemical reac- 
tions with electron transfer (ET) [17]. 

The Holstein model was originally proposed to charac- 
terize the vibration-assisted ET in one-electron molecu- 
lar crystal [17]. Here, in order to describe the chemical 
reaction of spin dependence as well as direction depen- 
dence, the Holstein model is generalized to incorporate 
the degrees of freedom of spin to make electrons natu- 
rally interact with a magnetic field. Additionally, due to 
the presence of the nuclear-spin environments surround- 
ing the electron spins, there would be coherent transi- 
tion between the singlet and triplet states of the two 
electron spins. In contrast to the previous investigation 
using anisotropic hyperfine coupling [11], the hyperfine 
interaction between the electron spin and its nuclear en- 
vironment is isotropic in our model. Based on this gener- 
alized model, we calculate the chemical reaction rates of 
the singlet and triplet states of the electron spins. Here, 
the chemical reaction rate is determined by the transi- 
tion rate of one electron in a localized molecular orbit 
to another at a distance. It is discovered that the re- 
action rate of the triplet state sensitively responses to 
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Figure 1: (color online). Schematic diagram for the gener- 
alized Holstein model. It is composed of a radical pair with 
nuclear spins. Two electrons are initially prepared in a cor- 
related state, i.e., the singlet state or triplet state, which can 
be interconverted by the hyperfine interaction in combination 
with an external magnetic field. The chemical reaction occurs 
once both electrons are in the same site due to the tunneling 
effect . 

the variation of the direction of the magnetic field with 
respect to the polarization of two electron spins. On 
the contrary, the chemical reaction of the singlet state 
does not demonstrate such dependence on the direction 
of the magnetic field. The above results are attributed 
to the invariance of the singlet state under the rotation 
of the system around ?/-axis, while the triplet one will be 
changed along with the rotation according to irreducible 
tensor of SO {3) group. Therefore, our proposed model 
may serve as a microscopic origin for the chemical reac- 
tion in the avian compass. 

In the next section, we generalize the Holstein model to 
incorporate the electron spin degrees. In Sec. Ill, we con- 
sider a general case with an external magnetic field and 
nuclear-spin environments. In Sec. IV, we study the dy- 
namic evolution of the radical pair and obtain the chem- 
ical reaction rates for different initial states. Finally, we 
summarize our main results in the Conclusion. Further- 
more, we show the detailed calculations for the chemical 
reaction probability, the chemical reaction rate and the 
transition probability from the triplet state to the singlet 
state in Appendix A and B respectively. 

II. GENERALIZED HOLSTEIN MODEL 

Many chemical reactions are accompanied by ET, 
where the electron is transfered from one site to an- 
other (Fig. 1). A very important but simple quantum- 



mechanical model for ET reactions is the molecular crys- 
tal model, which was originally developed by Holstein to 
describe so-called polaron motion in narrow-band con- 
ductors [18] and then understood as a microscopic model 
in describing simple chemical reactions [17]. 

The model Hamiltonian H = Hy + i^g + H^v is decom- 
posed into two parts, i.e., the ionic vibration part 

Hv = ^ ^(- — ^ + rrijuj'^jx]) + V^im2/xiX2,(l) 

j = l 3 3 

and the electron-phonon hybrid part 
He^Hev=^ej{xi,X2)c]^Cja-J ^(4^C2c + h.c.) , (2) 

where Xj is the displacement of the jth ion, uoj the har- 
monic vibration frequency with the reduced mass m^, / 
the coupling constant of the two molecules, and Cjoc{c^jo^) 
the fermionic annihilation (creation) operator of the elec- 
tron at molecule j with spin a. Since the orbital energy 
e{xi^X2) is linearized as 




we explicitly obtain the electronic Hamiltonian 

= ^^jc]a(^ja - J Y^{c\^C2a + h-C), (4) 

and the electron- vibration coupling 

Hey = — - AjC^-^CjgXj . (5) 

Here, the molecular orbital energy Sj is spin-independent. 
In the next section, we will consider a more general case 
with an external magnetic field. The tunneling inte- 
gral J is assumed to be independent of the displace- 
ment Xj and Aj = ^J^d'Sjldxj denotes the electron- 
vibration coupling, where the term dsj/dxi for i ^ j 
is neglected because the molecular orbital energy of the 
jth molecule changes negligibly when the displacement 
of the ith molecule varies. 

For simplicity, we assume two identical molecules, i.e., 
mi = m2 = m, uji =002^ and Ai = A2 = A. Choosing 
coordinates X = {xi -\- X2)/V^ and x = {xi — X2)/a/2, 
we decompose the Hamiltonian H = Hr + He into two 
decoupled parts, i.e., the one for the motion of the center 
of mass 

with Q = \/^\ + /, and the other for the relative motion 
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where we have introduced the bosonic operators 

^ = V^i ~ f effective frequency. 

Next, we make the Van Vleck transformation, also 
called the polaron transformation 0, HO], 

H = e'^He-'' (9) 

for the above generalized Holstein model, where 

s = -iY.<i>,{h'^ -h)c]^c,c. (10) 
is the transformation kernel and 

Thus we can formally decouple the degrees of freedom of 
electron and vibration, obtaining = Hr^^ + Hr^^ , with 

ja \ja J 

and 

F(i) = -J^(cLc2«e2*('''-'')+cLci«e-2*(''^-'')). (13) 

a 

This molecular crystal Hamiltonian ([12]) describes the 
ET process for a two-local-orbit system. Here, we gener- 
alize the Holstein model by taking into consideration the 
degrees of freedom of the electron spins. Up to now, the 
above generalization seems to be trivial, since we could 
totally separate the spin and orbital degrees of freedom. 
However, when a local external magnetic field is applied 
to the radical pair to form asymmetric couplings to the 
two electron spins, the spin-orbit coupling is induced. 
In this case, a spin-dependent ET process takes place. 
These asymmetric couplings can also be implemented by 
coupling to their nuclear-spin environments. 

III. SPIN MOLECULAR CRYSTAL IN 
MAGNETIC FIELD AND NUCLEAR 
ENVIRONMENT 




Figure 2: (color online). The inclination angle between the 
external geomagnetic field and the z-direction is G, and the 
coordinate system is rotated around the ^/-axis with G to co- 
incide with the direction of the the magnetic field. The hyper- 
fine coupling induced interconversion between the singlet and 
triplet states is modulated by the direction of the magnetic 
field, i.e., G. 

section, on account of an external magnetic field and 
nuclear-spin environments, we investigate how a chem- 
ical reaction responses to its magnetic environment. 

Choosing the polarization direction of the spin state as 
the z-direction, we define the singlet state \s) and triplet 
state \t) of the electron spins as 

i5)~(itU2)-iiit^))=^(4A-cI;4t)io), (14) 

and 

it)=^(im^)+ii?t^))=^(cVt^+4;4,)|o), (15) 

respectively with |0) being the vacuum state. 

In a simple case where the hyperfine couplings are 
isotropic, the Hamiltonian, which describes the inter- 
action between the electron spins and their asymmetric 
magnetic environments (the magnetic field plus the nu- 
clear spins), reads 



In the previous section, we described the generalized 
Holstein model, with spin degree of freedom. In this 
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Hs = - ^(/iB^o • Sj + gjij • Sj), (16) 
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where Bq = 5o(sin 6, 0, cos 9) is the external geomag- 
netic field with the inclination angle O, /i^ the Bohr 
magneton, Sj = (Sj^Sj^Sj) the Pauli operators for jth 

electron spin, Ij = (/J, /J, /J) the Pauli operators for jth 
nuclear spin, and gj the hyperfine coupling constant be- 
tween the jth electron spin and its environmental nuclear 
spin. 

Combining the relative vibration and spin Hamiltoni- 
ans, the total Hamiltonian for a spin-dependent ET re- 
action is obtained Htot = + Hs. After a rotation 
around ?/-axis with the angle 6 (Fig. 2), combined with 
the Van Vleck transformation defined in Eq. (|9j), we ob- 
tain Htot = Hr^Hs-, where the relative vibration Hamil- 
tonian Hr = Hr is the same as that given in Eqs. (|12til3p . 
but the Hamiltonian of the spin part is changed into 



Hs = — ^^^{I-IbBqSjz + Qjij • Sj). 



(17) 



Meanwhile, we make the same rotation and transforma- 
tion for the quantum states of the whole system as in 
Eq. Straightforwardly, after the combined transfor- 
mation, the singlet and triplet states read as \s) = \s) 
and 

\i) = cos e|i) + -i= sin einttl) - nut)) , (18) 

respectively. 

It is obvious that the singlet state is not B-dependent, 
while the triplet state is. In the next section, we will 
study the dependence of chemical reaction rate on the 
direction of the geomagnetic field. As shown as follows, it 
is this explicit dependence on the direction in the rotated 
triplet state that results in the variation of its chemical 
reaction rate along with the changes of the magnetic field 
direction. 



IV. MAGNETIC DIRECTION CONTROLLING 
CHEMICAL REACTION 

In this section, by means of the perturbation method, 
we analytically obtain the probability for one electron to 
transfer to the other local orbit to complete a chemical 
reaction. Assuming that at the initial time the vibration 
and nuclear spins are both in thermal equilibrium states 
and the electron spins are in the triplet state, the density 
matrix of the whole system p(0) = pv ® pn^ Pt includes 
three parts. The first part 



Pv 



^ CO 



exp 



m=0 



(19) 



denotes the relative vibration of the molecules, where 
Z = l/[l — exp{—tuj/kBT)], Ub the Boltzmann constant, 
and T the temperature of the environment. The second 
part is the density matrix of the nuclear spins. Since 
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Figure 3: (color online), (a) Chemical reaction probability vs 
time and G. {h) Chemical reaction probability for triplet state 
at a given time t = O.S/cj for effective temperature T = ^mK 
(red dashed line) and T = l^mK (blue solid line). 



PuBq <C ksT with pn the nuclear magneton, the nuclear 
spins are in the state 



^» = 7E|X")(X"I' 



(20) 



where |x?) = IIU?), Ix?) = lin?), Ix^) = Im2), and 
\xf) = |t?t2)- The last one 



Pt = \t) {t\ 



(21) 



describes the electrons. 

Starting from the above initial state, we calculate the 
total ET reaction probability of the triplet state (for the 
details please refer to Appendix A) 



oo 4 24 



^*(^) = ^ E EEe"'""'"^^-"^^' (22) 

where 

Pjmnp ('^) 



m,n=0 j = l p=l 



1 - e'' 



(23) 
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and Cjmq = ('0mg|?)|x?)l^) is the expanding coefficient. 



And the chemical reaction rate is determined by the reac- 
tion probabiUty per unite time in the long-time limit [2l| 



(a) 



d 

h = ^ lim Pt(r), 



(24) 



where the explicit expression which displays direction de- 
pendence is given in Appendix A. 

To show the above results in an intuitive way, we turn 
to numerical examples. We take the orbital energy dif- 
ference A = Si — €2 = O.OleV and the relative vibration 
frequency uj = lO^Hz. We assume the tunneling integral 
J = 0.01 A and (j) = 0.2. The magnitude of the geo- 
magnetic field is Bq = 50/iT and the hyperfine coupling 
constant is = ^2 = 10~^eV. In the following calcu- 
lations, we need to take a cut of the phonon occupation 
number, defined by an effective temperature of the envi- 
ronment. We find that the chemical reaction probability 
is not sensitive to the temperature of the environment. 
For these parameters we numerically calculate the tran- 
sition probability for initial triplet state from Eq. ([22]) as 
shown in Fig. 3. 

Obviously, the ET probability displays its dependence 
on the angle B. At a given time, the ET probability falls 
to its minimum value when the magnetic field is perpen- 
dicular to z direction, while it reaches its maximum when 
the direction of the external field is parallel to z axis. 
Besides, the probability is symmetrical about 6 = 7r/2. 
When we come to the initial singlet state, there is no such 
dependence on the angle O. This is a reasonable result 
for isotropic hyperfine coupling case since both the singlet 
state and the transformed Hamiltonian do not explicitly 
depend on the angle. 

On the other hand, the chemical reaction rate for 
triplet state would vanish once there were no interac- 
tion between the electron spins and their nuclear envi- 
ronments, i.e., gj = 0, while the ET reaction happens 
when the electron spins are in the singlet state. This can 
be seen from the fact that = but Hi^^^\s) ^ 0. 

The coupling of electron spins to the nuclear-spin envi- 
ronments can induce the transition from the triplet to the 
singlet states to make the ET happen. Then the spatial 
ET leads to the chemical reaction. In order to illustrate 
the mechanism of magnetic-direction-controlling chemi- 
cal reaction more clearly, we study the dynamic evolution 
of the radical pair. Using the same parameters as those in 
Fig. 3, we plot the transition probability from the triplet 
state to the singlet state 



Pt- 



-iHt/h~ 



me 



iHt/h 



(25) 



in Fig. 4. Due to the presence of the nuclear-spin en- 
vironments, there is periodical interconversion between 
the singlet and triplet states. Besides, the amplitude of 
this periodical oscillation is adjusted by the direction of 
the magnetic field. For one thing, the singlet-triplet in- 
terconversion is induced by the hyperfine couplings, not 
the uniform magnetic field, but it is modulated by the 
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Figure 4: (color online). Transition probability from the 
triplet state to singlet state, (a) Probability vs time for differ- 
ent angles, (b) The maximum probability varies with angle. 
The inset displays that the maximum transition probability 
also changes with the magnitude of the magnetic field. 



geomagnetic field. Our conjecture, that the triplet state 
is converted to the singlet state to complete the ET reac- 
tion, is conffimed by the same line shape of the maximum 
conversion probability as that of the chemical reaction 
probability of the triplet state. Besides, as shown in the 
inset of Fig. 4(b), with the increasing of the magnitude of 
the magnetic field, the maximum of the transition prob- 
ability increase at ffist for a given angle O = 0.1 tt and 
finally decreases into a stable value. 



V. CONCLUSION 

On account of electron spin degrees of freedom, we 
study the effect of the direction of the magnetic field on 
the chemical reaction by generalizing the Holstein model. 
By means of the perturbation approach, we obtain the 
ET reaction probability and chemical reaction rate of the 
singlet and triplet states. The chemical reaction rate of 
the triplet state displays its sensitive dependence on the 
direction of the magnetic field in contrast to the counter- 
part of the singlet state. We demonstrate that the triplet 
state indirectly participates in the chemical reaction. It 
must be converted to the singlet state by the hyperfine 
coupling between electrons and nuclear spins to take part 
in the ET reaction. We emphasize that the hyperfine cou- 
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plings are isotropic in our model which are different from 
the anisotropic ones in the previous study [11]. With 
the above comprehensive consideration, it could be con- 
cluded that our model may serve as a possible micro- 
scopic origin for the avian compass. 
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Appendix A: Chemical reaction rate 

In this appendix, we calculate the ET reaction proba- 
bility to the first order. 

First of all, the Hamiltonian of the electron spin-j part 
Hs^ is diagonalized with eigenstates 



e'^')=-o^im-\i)--^im^i)^ (A2) 



e«) = sin I \\.^) |t,"> + cos | |t,^) \V;) , (A3) 



Eigenstate 


Eigenvalue 












S,3 = eW+ef> 














1^7) = |e<^>) ® |ef > 


^,, = e^^) + ef) 


1^8) = |e«) ® |ef > 




l^g) = |e«) ® le*"') 


^,, = ei^) + e(^^ 


Ifio) = 14^*) <» |e^^*> 




\ipii) = |e^^*> |ef*> 




\ip12) = \ei^'') (» |ei^*> 




l^is) = \e[^^) <Si |ef *) 






i^sl4 — 64 +62 


IVis) = |ei^*> ® kf *> 


^sl5 — 64 +63 


Vi6> = e^^*> ® ef *> 




Iv'it) = nui) ® 14?;?) 


Esl7 = 


iv'is) = \mt) (s m^) 


^.18 = 


ivi9) = nut) ® It?;?) 


Esl9 — 


i¥'2o> = nut) ® itrt?) 


^s20 — 


b2i> = mm ® i;?;?) 


Es21 = 


l¥'22) = Iti^i) ® l^rt?) 


Es21 = 


l¥'23> = Iti^i) ® It?;?) 


^.22 = 


Iv'24) = Iti^i) ® itrt?> 


Es24 — 



Table I: All 24 eigenstates and eigenvalues of ifs, where the 
two electrons are located in two distant orbits respectively for 
the first 16 eigenstates, while for the other 8 eigenstates, both 
electrons are in the same site. 



e(^-)) = |t,«)|t,"), (A4) 

and the corresponding eigenvalues 

e^^^ = HBBo-gj, (A5) 
ei^^ = gj + yi^B^T^^ (A6) 

e^^ = g, - yi^B^T^^ (A7) 

ei^^ = -fiBBo-gj. (A8) 

Here, the mixing angles is defined as 

Straightforwardly, we explicitly calculate the eigenstates 
and eigenvalues of the total Hamiltonian for the spins 
Hs = H^s\ listed in Table I. 

The total Hamiltonian H = H^^'>^H^'^'> is split into two 
parts, i.e., H^^^ = ^ Hs governing the free dynamic 
evolution and H^^'^ = Hr^^ for the ET process. |V^mg) = 
\LPq) (8) \m) are the eigenstates of H^^^ with eigenvalues 



Erriq = TuhjO + Egq. Here, m (m = 0, 00) denotes the 
relative vibrational quantum number. A given initial state 
is expanded as 

24 

and then the wave function at time t is given by 

00 24 

= EE'^'^^We^pi-^^^^/^iK^)' (All) 

n=0 p=l 

where Cnp{t) are the coefficients determined by the 
Schrodinger equation 

Here, 

c.„p,„.,. = [<) - E'^lyn, (A13) 
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Coefficient 



Explicit expression 



+ sm 



,(0) 



+ cos^ ^sin^M(4%-<\o) 



- cos ^0 (^£,„_2l - ^™;9 



I + sin 



,6 



1 sin^ e[sin2 01 sin^ %5 (4°>3 - + sin^ 9^ cos* ^5 (^E', 



-Rl,mn,23 



n,23 ^m,7 



sin2e[cos^ ^sin^ 



+ sin^ 



r(o) 



„ L-v^^ 2 ^^^^ 2 ^ V ^n,21 ^m,2 / i ^ 

+ cos^ % cos* f<? «L-<\o)+ cos 



^m,3 /J 



R 



cos 



'2, mn, 22 



+ cos^ ^ sin' 



R2,nm,23 



+ sin^ ^1 sin^ 



i cos^ e[sin2 (?i sin^ 025 i - S^^^g ) + sin^ di sin^ $26 { 



(0) 
23 



p(0) 



AO) 



sin^ e[sin^ ^ cos^ «L - ) + cos' ^ sin^ - 
sin^ e[sin^ f cos^ 1 - ) + sin^ ^ cos^ f J (4'|>: 



R2,nn 



Rs,^ 



^ cos' e[sin2 ei sin' 626 [e^^^^ - E^°\) + sin^ Oi sin^ 628 (4°22 " 
+ sin^ e, sin' 626 «L - <\o) + sin^ 0: sin^ g^.? «L - <\i)] 



,mn, 22 



(0) 
m,7 



cos^ e[<5 (4% - <\0 + COS* f sin* f 5 (i?i°>3 " i^i^^e) + -s* f cos* f <5 (4% - 
+ sin* % sin* ^5 «>3 - + sin^ ^ cos* f ,5 (4% - E^, 



sin' e[sin' f cos^ (^"L - ^'m) + sin^ % cos^ f ^ - gj"',,) 



^3,1 



R4 



^[sin^ 01 cos* f 5 (<L - + sin^ 0^ sin* f 5 (<L - eI^],) 

sin' e, cos* f ^ «L - gi°\o) + sin^ 6, cos' ^& - 



R4,mn,22 



+ 



[cos^ ^ sin^ (<L - E^) + cos^ ^ sin^ 626 (<L - 
+ sin^ f sin^ fe^ «L - eZo) + ^in^ ^ sin^ 0,6 «L - ^i°!n 



Ra 



(0) 
m,7 



cos^ e[cos* (<L - <y + Bin* ^5 - eZi 

+ cos* f J «L - El^,^ + sin^ f J «L - i^i°\5) 



R4,mn,24 



Table II: Coefficients for the chemical reaction rate of the triplet state. 



and 



{m^q)^ we have 



(AM) 



To the first-order approximation, as Cn>p'{t) in the 
right-hand side of Eq. (|A12p is approximated as 
Cmp'(0)^n',m7 it is Straightforward to obtain 



Cnp(r)=cnp(0)-^c^,(0)^(2 



^ fiLUnp,mq 

(A15) 

The square of its norm gives the probability Pnp of finding 
the system in the state IV'np ) at time t. Thus, for {n,p) 7^ 



Pnpir) 



np,mq 



(AI6) 



The total ET reaction probability is P{r) = Pnp{r). 
And the chemical reaction rate is determined by the reac- 
tion probability per unit time in the long-time limit (2l|, 
i.e., 



d 

k = — lim P(r). 

OT r^oo 



(A17) 



In our case, the system is initially in the state p(0) 
Pv ^ Pn ^ Pti i-e., the relative vibration part 



1 



(AI8) 



m=0 
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Coefficient 


Explicit expression 


Dm,l,l 


1 [- sin ^ sin ^ cos Ge-^^-'^*/^ - cos ^ cos ^ cos Ge-^^-'^^/^ 
+ sin ^ sin ^ cos Ge-^^-'^*/^ + cos % cos % cos Ge-^^-'^'/^] 


-L^m, 1,2 


^ [- sin' ^ sin (92 sin Ge-^^-'^*/^ + sin' ^ sin 6*2 sin Ge-^^-'^*/^ 
- cos' ^ sin 6*2 sinGe"^^-^'io*/^ + cos' ^ sin 6*2 sin Ge"^^-'ii*/^] 


Dm,l,3 


^ [sin 6'i sin' ^ sin Ge"^^"^'6*/^ + sin Oi cos' ^ sin Ge"^^-'^*/^ 
- sin Oi sin' % sin Ge"^^-^'io*/^ - sin 6'i cos' % sin Ge"^^-'ii*/^] 


Dm,l,4 





Dm,2,l 


^ [- sin ^ cos f sin Ge-'^-'^*/^ + cos ^ sin ^ sin Ge-*^-'^*/^] 




1 [- cos Ge-^^-'4*/" + sin' ^ cos' f cos Ge-^^-'^*/" + sin' ^ sin' f cos Ge"^^-'^*/" 
+ cos' ^ cos' % cos Ge-^^-^'io*/^ + cos' ^ sin' % cos Ge-^^-^'^^*/^] 


Dm,2,3 


1 [- sin Oi sin 6*2 cos Ge"^^-'^*/^ + sin Oi sin 6*2 cos Ge"^^-'^*/^ 
+ sin 6'i sin O2 cos Ge-^^-'^o*/^ - sin Oi sin 6'2 cos Ge-^^-'i^*/^] 


Dm,2,4 


1 [cos ^ sin ^ sin Ge-^^-'^*/^ - sin ^ cos ^ sin Ge-^^-'^^t/zi] 


Dm,3,l 


^[sin ^ cos ^ sinGe-^^-'^*/^ - cos ^ sin ^ sin Ge-^^-'^*/^] 




1 [sin 6'i sin O2 cos Ge"^^-'^*/^ - sin Oi sin 6*2 cos Ge"^^-'^'/^ 
- sin Oi sin 6*2 cos Ge-^^'-'io*/^ + sin 6'i sin O2 cos Ge"^^-'ii*/^] 


Dm,3,3 


^ [- cos' f sin' f cos Ge-^^-'^*/^ - cos' ^ cos' f cos Ge"^^-'^*/^ 
- sin' f sin' f cos Ge-^^-'io^^^^o)'/^ - sin' ^ cos' f cos Ge-^^-'^^'/^ + cos Ge-^^-'^^V^] 


Dm,3,4 


1 [_ cos ^ sin ^ sin Ge-^^-'i4(v'i4)V^ + sin ^ cos ^ sin Ge-^^-'i^^'^i^)*/^] 


Dm,4,l 





Dm,4,2 


^ [sin 6'i cos' ^ sin Ge"^^-'6*/^ + sin Oi sin' ^ sin Ge"^^-'^*/^ 
- sin 6'i cos' ^ sin Ge-^^-'^o*/^ - sin Oi sin' ^ sin Ge-^^-'i^*/^] 


Dm,4,3 


I [- cos' ^ sin O2 sin Ge-^^-'^*/^ + cos' ^ sin 6*2 sin Ge-^^-'^*/^ 
- sin' % sin (92 sin Ge-^^-'^o*/^ + sin' ^ sin O2 sin Ge-^^-'^^*/^] 


Dm,4,4 


^ [- cos ^ COS f COS Ge-^^-'8*/^ - sin f sin f cos Ge-^^-'^^t/^ 
+ cos f cos f cos Ge-^^-'i4*/^ + gin ^ sin ^ cos Ge-^^-'^^*/^] 



Table III: Coefficients for the probability from the triplet to the singlet state at time t. 



with 



where 



Z = -exi){-huj/kBT)], 
the nuclear spin part 



Pn 



I 4 



and the electron spin part 

pt = \t) {t\ . 



(A19) 



(A20) 



(A21) 



After the Van Vleck transformation ([9]) and the ro- 
tation for the system, the initial state is transformed 
as p{0) = Pv ^ Pn ^ Pt, where py = py, pn = pn, and 
pt = |^)(^|- And we calculate the total chemical reaction 
probability as 

-J 00 4 24 

^*W = I7 E EEe"'""'"^™(^)' (A22) 



Pjmnp ('^j 



/ J ^jmq-LJ-np^nnq 



cjmq = {^i^l\i)\x])\rn). 



,(A23) 
(A24) 



As a matter of fact, those energy-increasing terms with 
q = 17, 18, 19, 20 (corresponding to the final states with 
both electrons at the orbits of the donor) and the cross- 
product terms with uJnp,mq 

7^ of Eq. (|A23l) do not con- 
tribute much to the chemical reaction probability. There- 
fore, a Fermi's golden-rule-like chemical reaction rate is 
obtained as 



^* = S^ E EE^"'""'"iHe-^*('''-''^HPi?. 

m,n—0 j=l p=21 



m,n=0 j—1 p=l 



with the coefficients being listed in Table IL 



jmnp 

(A25) 
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Appendix B: Singlet and triplet states 
interconversion 

For the system initially in the state 



(Bl) 



the probability of the electrons converted to the singlet 
state at time t reads 



(B2) 



4Z 



m,n 



X |(s|(x?l(n|e-'^*/'*|m)|xPW (B3) 

The eigenfunction of R is approximated to the first order 
as 

oo 24 

\^^,) = |^L°)) + E5]eKn,9,p)|Vi°)) (B4) 

n—0 p=l 



with 



^{m,n,q,p) 



^(0) _ ^(0) ' 



(B5) 



while the eigen energy is obtained to the second order as 



^mq - ^mq + (q) _ (Q) 



(B6) 



As a result, the time evolution operator is approximated 
as 

oo 24 
k=0 q=l 

Neglecting the the second-order terms, we obtain the con- 
version probability as 



4Z ^ ^ 

m=0 j,k=l 
16 

q=l 
^ oo 4 



4Z 



-jSmhoj It-) |2 
I J-^mj k I 7 



(B8) 



m=0 j,k=l 

where the coefficients Dmjk are listed in Table III. 
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